Substrate limited electron dynamics in graphene 
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We study the effects of polarizable substrates such as Si02 and SiC on the carrier dynamics in 
graphene. We find that the quasiparticle spectrum acquires a finite broadening due to the long- 
range interaction with the polar modes at the interface between graphene and the substrate. This 
mechanism results in a density dependent electrical resistivity, that exhibits a sharp increase around 
room temperature, where it can become the dominant limiting factor of electron transport. The 
effects are weaker in doped bilayer graphene, due to the more conventional parabolic band dispersion. 

PACS numbers: 71.20.Tx 72.10.Di 73.63.-b 



Introduction. Remote phonon scattering is a 
known limiting factor of the electron mobility in two- 
dimensional artificial structures of technological interest, 
such as Si MOSFETs [l|, i, S, i]. Due to the polar na- 
ture of the gate dielectrics used in such devices, the car- 
riers in the conducting channel couple electrostatically to 
the long-range polarization field created at the conduc- 
tor/dielectric interface. The interaction occurs primarily 
with optical phonon modes of high frequency (50 — 200 
meV), which do not affect transport at low temperatures. 
However, it results in a sizable degradation of the mobil- 
ity at room temperature (i.e. the temperature of device 
operation), and becomes the dominant scattering mech- 
anism in devices with high-K gate dielectrics. 

The interaction with the phonons of the gate dielectric 
has been recently investigated in another class of elec- 
tronic devices, where the Si inversion layer is replaced 
by a crystalline organic semiconductor 0, @] . In organic 
FETs, due to the narrow-band nature of the active mate- 
rial (typical bandwidths are less than 0.5 eV), the energy 
scale associated to the remote phonon interaction can be- 
come comparable with the bandwidth, resulting in carrier 
self-localization. In this case the effect is more dramatic 
than in Si MOSFETs: the mobility is not only strongly 
reduced by the interactions, but it also acquires an expo- 
nentially activated temperature dependence, due to the 
hopping motion of the localized carriers. 

Graphene -an atomically thin sheet of carbon atoms- 
is another system where the long-range coupling to the 
polar modes of the substrate can have a sizable influence 
on carrier dynamics. Compared to the two-dimensional 
electron gas formed in Si inversion layers, the effect here 
is enhanced due to the poor screening properties of the 
quasiparticles close to the nodal points, which behave 
as Dirac fermions. Although most experimental stud- 
ies of transport in graphene up to now have focused on 
the low temperature regime [3, II], its possible use as 
a material basis for future electronic devices calls for a 
deeper understanding of the room temperature behavior 



In this work we calculate the effects of remote phonon 
scattering on the dynamics of electrons in graphene [l^ . 
Unlike other scattering mechanisms such as disorder 
[isl . [3, H, 18, [l^, or acoustic phonon scattering [20| . 
all the microscopic parameters of the remote electron- 
phonon interaction are known a priori from indepen- 
dent measurements, so that in principle its effects can 
be estimated quantitatively. We show that in epitax- 
ial graphene grown on SiC, the polar phonon scatter- 
ing has only a weak effect on the electron mobility, due 
to the weak polarizability of the substrate and the rel- 
atively high phonon frequencies associated to the hard 
Si-C bonds. The effect is much larger at graphene/Si02 
interfaces where, at room temperature and at the typical 
densities n ^ IO^^cto"^ attained in current experiments, 
it gives rise to resistivities of the order of lOOfi. This sug- 
gests that, as in traditional Si based electronic devices, 
remote scattering with the substrate phonons can indeed 
constitute an important limiting factor of the mobility in 
future graphene devices. 

The model. The polar phonons of the substrate in- 
duce an electrostatic potential which couples to the elec- 
tron charge. This potential is smooth over lengths com- 
parable to the lattice constant, due to the finite distance 
between the substrate and the graphene layer (see be- 
low) , which allows us to neglect intervalley scattering. In 
this case the electron-phonon interaction term is given 
by Hi — ^qMqPq{hq + fol^), where pq is the electron 
density operator and the operators 6 J , bq describe the 
phonon displacements. The interaction matrix element 
can be written as ^ ge^'^iyiqa) UMM^, where z 
is the distance between the graphene layer and the sub- 
strate and a = 1.42 A is the lattice spacing. Neglecting 
the dielectric response of the atomically thin graphene 
layer, the dimensionless coupling parameter g is given by 



5 = 27r/3 



hvp/ahvp 



(1) 



Here /? — {cs — f-oo)/{f-s + l)/(eoo + 1) is a combination 
of the known dielectric constants of the substrate, vp = 



10^ cm/s the Fermi velocity and utg the frequency of the 
interface phonon mode under consideration. We consider 
separately the case of epitaxial graphene grown on SiC, 
and that of graphene flakes deposited on Si02 substrates. 

In 6H-SiC there is a single surface phonon mode at 
LUs — 116meV (^^l- The dielectric constants of the bulk 
material are = 9.7, eoo = 6.5, which yields f3 = 0.040 
and g — 1.7 • 10"^. In the case where several surface 
modes are present, as in Si02, each mode will couple 
through a partial factor /3i weighted by its oscillator 
strength (the sum of all Pi's still given by (3 defined 
above). For example, in crystaUine Si02 (cs = 3.9, 
Eoo = 2.4), there are two dominant surface modes at 
LUsi = 59 meV and uJs2 = 155 meV, with (3i = 0.025 
and (32 = 0.062 respectively (/3 = 0.087) [i], correspond- 
ing to 3i = 5.4-10"'^ and 52 = 3.5-10~^. These values are 
enhanced by roughly 50% in common Si02 glass, where 
Eoo ~ 2.1. For completeness, we also give the values 
for the high-K dielectric IIf02 that has been used in re- 
cent experiments on locally gated graphene [l^: = 22, 
= 5, /3 = 0.12 and ~ MmeV , so that g = 4.2- 10"^ 

Scattering rate. The quasiparticle scattering rate 
arising from the remote electron-phonon interaction de- 
scribed above can be calculated as 

TH = |^^p,qA/|{(5(ep+q- |w + w,|)K-|-n+] 
q 

+ 5{ep+^-\uj - ujs\)[nB + 1 - Up]] , (2) 

where the band dispersion is a; = ±ep ~ ±vfP, ub is the 
Bose distribution for phonons of energy ojs and , Up 
are the Fermi functions for electrons at oj ± oj^. The fac- 
tor Ap q = [1 + scos((/)p+q — 0p)]/2 is the spinor overlap 
for intraband (s = 1) and interband (s = —1) scattering 
respectively, where (/)p defines the direction of p. To a 
first approximation, dynamic screening from the conduc- 
tion electrons can be accounted for by replacing l/g — > 
1 / (q + Qs) in the interaction matrix element M^, where 
Qs = ^e^kp /hvp is the Thomas- Fermi screening length 
[2J]. We set the graphene/substrate distance to the typ- 
ical value z = 4A. Eq. 12]) can be directly generalized to 
bilayer graphene, by considering a parabolic band disper- 
sion tp = Vpp^ /t±, with t± ~ 0.35ey the interlayer hop- 
ping parameter. In that case the screening wavevector 
becomes independent of band filling, = At±/hvF, and 
the spinor overlap changes to [1 -I- s cos 2(0p+q — (f>p)]/2. 
We also take a slightly larger separation z ~ 6A, cor- 
responding to the average distance of the bilayer to the 
interface, although in principle the electron-phonon cou- 
pling differentiates between the two graphene sheets. 

In Fig[T] we show plots of the scattering rate F cor- 
responding to single-layer and bilayer graphene on the 
two different substrates SiC and Si02, for both undoped 
{Ep ~ 0) and heavily doped {Ep = 0.6eV) graphene. 
In all cases the scattering rate is suppressed in a win- 
dow ±LOs around the Fermi energy (it actually vanishes 
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FIG. 1: Quasiparticle's scattering rate due to remote phonon 
interaction as function of energy for: single-layer (top) and 
bilayer (bottom) graphene, with the parameters appropriate 
to an SiC (left) and a Si02 (right) substrate respectively. 
The full and dotted lines correspond respectively to undoped, 
Ep = 0, and heavily doped, Ep = 0.6eV cases. Notice the 
different axis scale in single-layer (top) and bilayer graphene 
(bottom). 

at T = 0, because in that case phonon exchange is for- 
bidden in all this range) . The absolute value of F at the 
chemical potential, that will enter in the calculation of 
the d.c. conductivity, is proportional to both the cou- 
pling strength g and the phonon thermal population, i.e. 
F cx 5e~"=/^. Note that the anomalously large scatter- 
ing rates obtained for undoped single-layer graphene (full 
lines in the upper panels of Fig. 1) compared to all the 
other cases are related to the poor screening properties 
of Dirac electrons close to the neutrality point. 

The present results for the quasiparticle lifetimes from 
remote phonon scattering can be compared with analo- 
gous estimates for the interaction with the intrinsic in- 
plane vibrations of the graphene sheet [25[ . The effect of 
the substrate appears to be slightly weaker than the in- 
trinsic electron-phonon scattering in the case of epitaxial 
graphene grown on SiC, while the two mechanisms are 
of comparable magnitude at graphene/Si02 interfaces. 
Nevertheless, we will show below that due to the smaller 
oscillation frequencies involved, remote phonons have a 
much stronger influence on electronic transport than the 
intrinsic phonons of graphene. 

Conductivity. We calculate the electrical conductiv- 
ity through the Boltzmann equation. Taking into ac- 
count spin and valley degeneracy, the d.c. conductivity 
of single-layer graphene can be expressed as: 

a=^|d.MFr/(c.)(-^). (3) 

Here Ttr{i^) is the transport scattering rate, defined 
as in Eq.([2|), with the additional angular factor [1 — 
s cos((/)p+q — (j)p)] in the integrand which favors large an- 
gle scattering events [l^. Note that this approximation 
is strictly valid only in the elastic limit ivq —>■ (where 
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Ttr ~ r/2) while it slightly underestimates the scattering 
rate in the inelastic regime co ^ loq- In the case of bilayer 
graphene, the integral in Eq. ([3]) has an additional factor 
2uj/t± arising from the parabolic band dispersion. 

We now present results for graphene on Si02 sub- 
strates, where most experimental transport measure- 
ments have been performed (the predicted effect at room 
temperature is one order of magnitude smaller in the case 
of SiC). From Eq. ([3]), since the derivative of the Fermi 
function selects contributions close to the Fermi energy, 
we see that the conductivity is inversely proportional to 
the scattering rate at ~ Ep, which in turn scales as 
g^-ujs/T j-ggg above). Since the coupling constant g in 
Eq. ([1]) is directly proportional to the phonon frequency 
LUs, we see that the scattering rate at a given temperature 
T is a non-monotonic function of the phonon frequency: 
it is maximum at lUs = T, and decreases exponentially 
for larger phonon frequencies. As a result, the phonons 
that mostly affect transport will be the ones with the 
lowest frequencies (closest to T). In particular, in the 
case of a Si02 substrate, the dominant phonon is the 
one at ujsi = 59meV, even though it has a weaker cou- 
pling strength (due to the exponential temperature de- 
pendence, scattering from the phonon at ujs2 — 155meV 
is mostly inactive at room temperature). For the same 
reason, in-plane optical phonons whose characteristic en- 
ergy scale is ^ 200meV should give a negligible contribu- 
tion to the resistivity at room temperature. 

Results. In Fig. [2K we show the result of p vs T 
at Ep = lOOmeV (corresponding to a density n = 
0.7 • lO^^cm"^) taking a finite distance z = iA between 
the substrate and the graphene layer and a Thomas- 
Fermi screening from the graphene electrons as described 
above [i^]. As anticipated from the preceding discussion, 
the resistivity exhibits a characteristic exponential be- 
havior. Note that unlike the short-range scattering from 
in-plane phonons, the present result has a sizable density 
dependence due to the long-range nature of the inter- 
action. This is illustrated in Figllb, where we plot the 
dependence of a on Ep at fixed T = 300, 200, 150K (from 
bottom to top). In the explored density range, the con- 
ductivity has an overall linear dependence on Ep away 
from the Dirac point, which is intermediate between the 
cases of short range scatterers (cr ~ const) and charged 
impurities {a cx _E^)[l5l|. This behavior can be under- 
stood from standard dimensional arguments, observing 
that the scattering rate of Eq.([2]) scales as M^^^^ oc kp^ 
times the density of states cx kp, which yields T ~ const. 
The result a cx \Ep\ then follows directly from Eq. 

Superimposed to the smooth density dependence an- 
alyzed above, a slight "kink" appears in the electrical 
characteristic at Ep ~ cj^i that testifies the interaction 
with the sharp interface phonon mode. The anomaly be- 
comes more pronounced at low temperatures owing to the 
sharp Fermi function in the Kubo formula. There, how- 
ever, the absolute value of the scattering rate becomes 
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FIG. 2: a) Temperature dependence of the resistivity for 
Ep = 100 meV (the dashed line corresponds to p = poe~"°^''"^ 
with po ~ 50051). b) Conductivity vs. chemical potential for 
T = 150, 200 and 300K (from top to bottom). The inset of 
a) shows the comparison with bilayer graphene (dotted line) 
at T = SOOK. 



extremely small, and it is likely to be completely hidden 
by the presence of other sources of scattering. 

In the inset of Fig. [^b we compare the conductiv- 
ity for single-layer and bilayer graphene at T = 300K. 
The different density of states in bilayer graphene leads 
to a stronger energy dependence. Accordingly, bilayer 
graphene should be a better conductor at large bias. 
With the present choice of parameters, the crossing oc- 
curs at Ep ^ O.leV^. 

Finite distance cutoff. The simple scaling argument 
presented above, leading to a linear dependence a oc \Ep\ 
for single layer graphene, breaks down when 2k pz ~ 1. 
At larger densities, the short range cut-off associated to 
the distance z between the carriers and the interface must 
be taken into account explicitly. As a result, extra powers 
of kpz arise in the scattering rate, due to the suppres- 
sion of the angular factors [1 ± cos((/)p+q — (j)p)] at large 
scattering angles, as well as the reduced available density 
of states. All together, these factors eventually lead to 
a (7 oc Ep dependence (F oc |i?i?p). In Fig. 3 we show 
plots of cr vs Ep a.t T — 300K for increasing values of z. 
While for z — 4 A the behavior is approximately linear 
in the whole energy range, a clear upturn becomes visible 
at larger separations, reflecting the expected crossover at 
\Ep\z/vp ^ 1/2. Notably, this phenomenon can lead to 
an apparent behavior a oc Ep oc n Ts'l, which is best 
seen in the inset of Fig. 3. 
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FIG. 3: Same as Fig 2.b, for different values of the distance 
z. From bottom to top, z = 4, 12 and 20 A. A crossover to a 
super-linear behavior is apparent at large z and large Ef- In 
the inset, the same data are plotted as a function of oc n. 

Discussion. Although the transport properties of 
graphene are roughly temperature independent below 
T = lOOK, a sharp increase of the resistivity of the or- 
der of Ap ^ lOOri around room temperature has been 
reported in field-effect doped samples at n > lO^^cm"-^ 
0, [13, Ell, [S]- This increase is unlikely to be due 
to the effects of disorder, which are expected to depend 
smoothly, if at all, on the temperature. As was mentioned 
above, also the in-plane optical phonons of graphene can 
be ruled out due to their high oscillation frequencies. 
That in-plane scattering is not the main cause of the 
observed behavior is further supported by the observa- 
tion that the temperature dependent contribution to the 
resistivity is not independent of density [ill [isj , point- 
ing to a long-range scattering mechanism such as the one 
studied here. According to the present results, a density 
and temperature dependent contribution to the resistiv- 
ity is compatible with the mechanism of remote phonon 
scattering. 

Conclusions. We have estimated the quasiparticle 
scattering rates induced by the polar modes of Si02 and 
SiC substrates on graphene. Our results show that, while 
remote phonon scattering has negligible effects on trans- 
port at low temperatures, it becomes relevant at room 
temperature, where it could constitute the dominant lim- 
iting factor of the electron mobility in sufficiently clean 
samples, especially for Si02 substrates. According to the 
present scenario, the maximum conductivity that can be 
achieved in a given sample depends crucially on the di- 
electric properties of the gate insulator, so that the use 
of non-polar substrates should be favored to optimize 
the transport properties of graphene. We also find that 
the effects are reduced in doped bilayer graphene, which 
makes it more suitable for device applications at room 
temperature. On the opposite, larger effects can be ex- 
pected when more polarizable substrates are used as gate 
dielectrics. 
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